arXiV: astro-ph/0112412v 1 17 Dec 2001 


Faraday Rotation Effect of Intracluster Magnetic Field on Cosmic Microwave 

Background Polarization 

Masahiro Takada^’^, Hiroshi Ohno^, and Naoshi Sugiyama^ 
mtakada@th.nao.ac.jp; ono@resceu.s.u-tokyo.ac.jp; naoshi@th.nao.ac.jp 

ABSTRACT 

The observed magnetic field of microgauss strength in clusters of galaxies should induce the 
Faraday rotation effect on the linearly polarized cosmic microwave background (CMB) radiation 
when the CMB radiation propagates through a cluster at low redshift. The Faraday rotation arises 
from combined contributions of the magnetic field strength, the electron density, the cluster size, 
and the characteristic scale of the magnetic field along the line of sight. Employing the Press- 
Schechter prescription for the cluster abundance under the cold dark matter (CDM) scenario 
and a plausible isothermal /3-model for the gas distribution, we compute angular power spectra 
of the CMB polarization fields including the Faraday rotation mixing effect under the simple 
assumption of uniform magnetic field configuration across a cluster. As a result, we find that a 
parity-odd B-type polarization pattern is statistically generated on the observed sky, even when 
the primary polarization only contains the parity-even E-mode component, such as in the case 
of pure scalar perturbations. The generated E-type polarization has a peak with amplitude of 
^ 0.1 /rK(Eo/0.1 /rG)(t'o/10 GHz)“^ at angular scales of I « 1000 for the currently favored 
adiabatic ACDM model. This result also implies that, if the magnetic field has a 0.5 p,G strength 
and we observe at lower frequencies such as izq ^ 5 GHz, the secondary signal due to the Faraday 
rotation effect could be comparable to the magnitude of the primary polarization. The frequency 
dependence of the Faraday rotation can be then used to discriminate the effect from primary and 
other secondary signals on the CMB polarizations. Our results therefore offer a new empirical 
opportunity to measure or constrain the intracluster magnetic field in the average sense, combined 
with measurements of the intracluster gas distribution through the X-ray and SZ data. 


Subject headings: cosmology: theory - cosmic microwave background - cluster of galaxies - 
magnetic fields - polarization 


1. Introduction 


Various observations have revealed that clusters of galaxies are pervaded by the strong magnetic field 
of microgauss strength (e.g., see Kronberg 1994; Carilli, & Taylor 2001 for reviews and references therein). 
The multi-frequency Faraday-rotation measurements of polarized radio sources inside or behind a cluster 
have been used to estimate the magnetic field strength, combined with the X-ray data (Vallee, MacLeod, fc 
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Broten 1986; Kim, Tribble, & Kronberg 1991). Recently, Clarke, Kronberg, & Behringer (2001) have drawn 
a firm conclusion that an intracluster hot plasma is universally magnetized by 1 — 10/rG fields under the 
assumption of 10 — lOOkpc field coherent scales, using 16 normal low-z ROSAT cluster sample selected to be 
free of unusual strong radio halos, widespread cooling flow and recent merger activity. However, except for 


a few cases such as some background rotation measures per square degree in the Coma cluster (Kim et al 


1990|), it is generally difficult to measure the angular profile of the magnetic field strength inside a cluster 


because of the lack of the number of radio sources available per cluster. The observations of cluster-wide 
diffuse radio halos have also led to the evidence of the microgauss field, which is believed to be synchrotron 
radiation by relativistic electrons accelerated in the shock wave. Moreover, in this case hard X-ray emission 
could be produced by the same relativistic electron population through the inverse Compton scattering 
of cosmic microwave background (CMB) photons. The combined observation of hard X-ray emission and 
radio halos is therefore very attractive in the sense that it allows us to directly estimate the magnetic field 
strength without further restrictive assumptions for the coherent length and the electron distribution. The 
recent detections of hard X-ray emission have led to an independent estimation of ~ 0.1/xG fields even in 
the outer (1 Mpe) envelopes of clusters (Bagchi, Pislar, fc Neto 1998; Rephaeli, Gruber, fc Blanco 1999|; 


Fusco-Femiano et al. 1999). 


However, the origin of the intracluster magnetic fields still rests a mystery in cosmology. The following 
two scenarios have been usually investigated in the literature. One is based on the idea that the cluster 
field is related to the fields generated by dynamo mechanism in individual galaxies and subsequent wind-like 
activity transports and redistributes the magnetic fields in the intracluster medium (e.g. Kronberg, Lesch 


fc Hopp 1999 and references therein). The general prediction of galactic dynamo mechanism is that the 
galactic magnetic field could arise from an exponential amplification of a small seed field during a galactic 
lifetime and energy of the mean magnetic field could grow up to equipartition level with the turbulent energy 


of fluid (Parker 1979; Zel’dovich, Ruzmaikin, & Sokoloff 1983; Lesch, fc Chiba 1997). It is not still clear, 
however, that the galactic dynamo theory can explain the detection of microgauss magnetic fields in high-z 


damped Lya absorption systems, which are supposed to be protogalactic clouds (Welter, Perry, fc Kronberg 


1984; Wolfe, Lanzetta, fc Oren 1992). On the other hand, an alternative scenario is that the intracluster 


magnetic fields may grow via adiabatic compression of a primordial field frozen with motions of cosmic 
plasma, where the primordial field is assumed to be produced somehow during the initial stages of cosmic 


evolution (e.g., Rees 1987; Kronberg 1994; and also see Grasso, fc Rubinstein 200C for a recent review). 
Recently, using the cosmological, magneto-hydrodynamic simulations of galaxy clusters, Dolag, Bartelmann, 
fc Lesch (1999) have quantitatively shown that, if starting with the initial magnetic field of ^ 10“®G strength, 
the final intracluster field can be amplified by the gravitationally induced collapsing motions to the observed 
strength of ~ fj,G irrespective of uniform or chaotic initial field configurations motivated by scenarios of the 
primordial or galactic-wind induced initial seed fields on Mpe scales, respectively. It is undoubtedly clear 
that the knowledge of the intracluster magnetic field will lead to a more complete understanding of the 
physical conditions of the intracluster medium and of the possible dynamical role of magnetic fields (e.g.. 


Loeb, & Mao 1994; Vikhlinin, Markevitch, & Murray 2001). Furthermore, the magnetic field should play an 
essential role to the nonthermal processes such as the synchrotron and high energy radiations and possibly 


the cosmic ray production ( Loeb, fc Waxman 2000 ; Totani, fc Kitayama 2000 Waxman, fc Loeb 2000 ). 
Hence, it is strongly desirable to perform further theoretical and observational investigations on these issues 
in more detail. 

In this paper we investigate how the Faraday rotation effect due to the intracluster magnetic field causes 
a secondary effect on the CMB polarization fields. While Thomson scattering of temperature anisotropies 
on the last scattering surface generates linearly polarized radiation at the decoupling epoch ( Kosowsky 1996| ; 
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Kamionkowski, Kosowsky, fc Stebbins 1997a ,b; Zaldarriaga, fc Seljak 1997 ; Hu, fc White 1997^ b), the plane 
of linear polarization should be rotated to some extent when the CMB radiation propagates through the 
magnetized intracluster medium at a low redshift. The great advantage is that, since the CMB polarization 
field is a continuously varying field on the sky, we can measure an angular profile of the rotation measure in 
a cluster in principle, which cannot be achieved by any other means. Furthermore, thanks to the frequency 
dependence of the Faraday rotation, one will be able to separate this effect from the primary signal and other 
secondary signals such as that induced by the gravitational leasing effect due to the large-scale structure 


( Zaldarriaga, fc Seljak 1998 ). In this paper, based on the Press-Schechter theory for the cluster abundance 
under the cold dark matter (CDM) structure formation scenario and a plausible isothermal /3-model of the 
intracluster gas distribution, we compute angular power spectra of the CMB polarization fields including the 
Faraday rotation mixing effect caused by clusters at low redshifts. As for the unknown field configuration, 
we assume the uniform field of ~ 0.1/rG strength across a cluster that is consistent with the recent Faraday 
rotation measurements (Clarke et al. 2001). This model allows us to estimate the impact of this effect in the 
simplest way. Our study thus proposes a new empirical opportunity to measure or constrain the intracluster 
magnetic field in the average sense that properties of the magnetic field could be extracted through changes 
of the statistical quantities, CMB angular power spectra. 

So far previous works have focused mainly on investigations of effects of the primordial magnetic fields 
on the CMB temperature and polarization anisotropies (Kosowsky, & Loeb 1996; Aduins_e^_aL^199^ 


Scannapieco, & Ferreira 1997; Harari, Hayward, & Zaldarriaga 1997; Seshadri, & Subramanian 2001; Mack 


Kahniashvili, fc Kosowsky 2001 ; and see also Grasso, fc Rubinstein 200C for a review). Both the amplitudes 
of the cosmological primordial magnetic field and the mean baryonic density rapidly increase with redshift 
scaled as (1 -I- z)^ and (1 -I- z)^, respectively. The dominant contribution to the Faraday rotation effect is 
therefore imprinted before decoupling (z ^ 10^), unless the structure formation at low redshifts causes a 
significant amplification of the magnetic field. Kosowsky, & Loeb (1996) found that, if the current value of 
the primordial field is of the order of 10“®G on Mpc scales corresponding to the microgauss galactic field 
under the adiabatic compression, the effect on the CMB polarization is potentially measurable by satellite 


missions MAP and Planck Surveyor (see also Bcannapieco, fc Ferreira 1997; Harari, Hayward, fc Zaldarriaga 


1997). It is also shown that such a stochastic magnetic field in the early universe affects the temperature 


fluctuations through the induced metric vector perturbations and thus current measurements can put an 
upper limit of ;$ 10“^^G on the current field amplitude (Mack, Kahniashvili, & Kosowsky 2001). 


This paper is organized as follows. In Section g we first construct a model to describe the magnetized hot 
plasma in a cluster by assuming the uniform magnetic field configuration for simplicity. We then derive the 
angular power spectrum of the Faraday rotation angle based on the Press-Schechter description of the cluster 
formation and the /3-model of the gas distribution. In Section we present a formalism for calculating the 
angular power spectra of the CMB polarization fields including the Faraday rotation mixing effect. Section 
^ presents the results for the currently favored CDM models. In Section ^ we briefly present the discussions 
and summary. Unless stated explicitly, we assume the favored ACDM cosmological model with Hmo = 0.3, 
Hao = 0.7, = 0.05, h = 0.7, and cts = 1-0 as supported from observations of CMB anisotropies and 

large-scale structure (e.g. Netterfield et al. 2001), where Hmo, and Had are the present-day density 
parameters of non-relativistic matter, baryon and the cosmological constant, respectively, h is the Hubble 
parameter and cts denotes the rms mass fluctuations of a sphere of 8h“^Mpc radius. We also use the c = 1 
unit for the speed of light. 
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2. Model of Faraday Rotation due a Cluster Magnetic Field 
2.1. Faraday rotation of a magnetic field 


If monochromatic radiation of frequency v is passing through a plasma in the presence of a magnetic 
field B along the propagation direction 7 , its linear polarization vector will be rotated through the angle 
Atp (e.g., Rybicki, fc Lightman 1979): 


Ai^ = 


2'Krai(?v^ 


JdlueiB ■ 7 ) 


8.12 X + z) 


Aq 

1 cm 


dl 

kpc 


( ) ( 

V 1 cm“3/ 


B 


7 


1 fj.G 


( 1 ) 


where e and me denote the charge and mass of electron, respectively, Ue is the number density held of 
electron along the line of sight, and Aq is the observed wavelength. Note that the factor (1 + z)~'^ comes 
from the relation of A = Ao(l + 2 )”^ and the sign of A(p could be negative depending on the orientation 
of the magnetic held and the reversal scale of magnetic held. The equation above clearly shows that the 
magnitude of the rotation measure depends on the electron number density, the magnetic held strength 
and the characteristic scale length of the held, and is larger for longer observed wavelengths because of the 
dependence Aq. Thus, as a possible source of the Faraday rotation, in this paper we consider the magnetized 
hot plasma in clusters of galaxies, motivated by the fact that the deep gravitational potential well of dark 
matter makes the intracluster gas dense and highly ionized and many observations have revealed the existence 
of the strong magnetic held with ~ fiG strength(e.g., Clarke et al. 2001). 


2.2. Cluster model with a uniform magnetic field 

To calculate the Faraday rotation measure for the magnetized intracluster hot plasma, we have to 
assume how the magnetic held is distributed in a cluster along the line of sight. There are two possibilities 
considered; the uniform held conhguration and, more realistically, the tangled held, where the former held 
has constant strength and direction through the entire cluster while the latter held has reversal scales along 
the line of sight. For the Coma cluster, Kim et al. (1990) suggested ~ 10 kpc for the held reversal scale 
from the observed rotation measures of some sources close in angular position to the cluster, and Feretti et 
al. (1995) discovered the magnetic held structure with smaller scales down to ~ 1 kpc from the rotation 
measures of radio halo in the Coma, leading to the stronger held of ^ 10/xG on the scale. On the other 
hand, there is observational evidence of coherence of the rotation measures across large radio sources, which 
indicates that there is at least a held component smoothly distributed across cluster scales of ~ 100 kpc (e.g., 
Taylor, fc Perley 1993| ). Thus, the magnetic held structure is still unknown because of the lack of the number 
of sources available per cluster for measuring the magnetic held. The angular prohle of the rotation measure 
compiled from all sample of radio sources for 16 different normal clusters (Clarke et al. 2001) yields the held 
strengths between ~ 0.5 and ~ 3 fxG for a uniform held model, while ~ 5 fiG{l/10 kpc)“^/^(/i/0.75)^/^ for a 
simple tangled-cell model with a constant coherence length 1. In this paper, we consider a uniform magnetic 
held model for the sake of revealing the impact of the Faraday rotation effect on the CMB polarization in 
the simplest way. We then assume that the amplitude of magnetic held at the cluster formation epoch is 
universally constant, namely it does not depend on redshift; i3(zform) = Bq = constant. 

Next let us consider a model to describe the gas distribution in a cluster. Motivated by the X-ray 
observations, we adopt the spherical isothermal /3-model as for the gas density prohle of each cluster expressed 
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by 


p{r) = po[l + (r/rc)'^] , ( 2 ) 

where po is the central gas mass density and is the core radius. The observed values of j3, obtained from 


X-ray surface brightness profiles, range from 0.5 to 0.7 (e.g. see Jones, & Forman 1984). In the following 
discussion we assume (3 = 2/3 for simplicity. Note that Makino, Sasaki and Suto (1997) showed that the 


universal density profile of dark matter halo (Navarro, Frenk, & White 1997) can reproduce the /3-model 
for massive clusters under assumptions of the hydrostatic equilibrium and isothermality. In the spherical 
collapse model, the virial radius rvir for the halo with mass M is obtained from 


M = yPcr(^)A(z)r^ir, 


( 3 ) 


where A(z) is the overdensity of collapse and the fitting formula is given by Bryan & Norman (1998), and 
Pcr{z) is the critical matter density at redshift z defined by Pcr{z) = “HIq /+ z)^. 

Once the uniform magnetic field and the gas density profile are given, we can easily calculate the angular 
profile of the Faraday rotation angle for a cluster projected on the sky similarly as the procedure to calculate 


the profile of the thermal Sunyaev-Zel’dovich (SZ) effect ( Atrio-Barandela, fc Mficket 199E ; Komatsu, fc 
Kitayama 1999). Using the gas density profile (^) and equation (^, the angular profile can be calculated as 


A(p{e, 


e'^Ueorc 


27rm2c2j/2(i + 2)2 


{B cos x) 
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( 4 ) 


where p = r^irlrc, X is the relative angle between the direction of the magnetic field and the line of sight, 
rieo is the central number density of electron, and 9 is the angular separation between the line of sight and 
the cluster center. The angular core size 9c is related to via 9c = rc/d^^^z), where dA{z) is the angular 
diameter distance to the cluster at redshift z. Note that p, given as the ratio between the viral and the core 
radius, is not a parameter but a function of mass M and redshift z in a general case. 

We also have to relate the quantities Tc and rieo to the total mass M and redshift z. Assuming the full 
ionization for intracluster gas, the central electron number density is given by rieO = (Pgaso/wp)(l — Y/2), 
where Y is the helium mass fraction and we fix U = 0.24 in this paper. It is known that the formation of 


cluster core regions contains large uncertainties. We here employ the entropy driven model ( Kaiser 1991 ; 
Evrard, & Henry 1991) to describe the core radius, motivated by the fact that this model can explain 


observations such as the X-ray luminosity function better than the simplest self-similar model ( Kaiser 1986| ) 
does. In this model the core is supposed to be in a minimum entropy phase Smin = c„ ln(Teo/Pg//o)j where Tco 
denotes the central electron temperature, Cv is the specific heat capacity of the gas at constant volume, and 7 
is the ratio of specific heats at constant pressure and constant volume. If using M = ( 47 r/ 3 )pcr(z)A(z)r^ij. « 
47rr/rvir(pgaso//gas) for Tvir/rc > 1 , we can obtain 
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where /gas is the gas mass fraction for a cluster and we assume that /gas is equal to the cosmological mean; 
/gas = Mgas/M = Ubo/^^mo (White et al. 1993). The adiabatic index is fixed at 7 = 5/3 throughout this 


paper. Moreover, we assume the isothermality and that the central electron temperature Tgo is equal to the 
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virial temperature Tvir given by 


ksTyir = 77 


fimpGM 

Sr ■ ’ 

01 vir 


( 6 ) 


with /r = 0.59. The quantity 77 is a fudge factor of order unity that should be determined by the efficiency of 


the shock heating of the gas (Eke, Navarro, & Frenk 1998), and we fix 77 = 3/2. If the proportional constant 
to calculate Vc in equation (||) is determined by assuming rc{10^^h~^M q, z = 0) = 0.15/i“^Mpc, we can 
obtain the core radius for an arbitrary cluster with mass M and at redshift 2 ;. This cluster model then yields 
Tvirlrc = 13.6, Ueo = 5.0 X 10“^ cm“^, and Te = 6.5keV for M = Mq at z = 0 in the ACDM model. 

Using these typical values of physical quantities for a cluster with M = lO^^/i“^M 0 , equation i) tells us 
the magnitude of the rotation measure at the cluster center as 


A(77(z = 0,0 = 0) 


2.34 X 10"^ rad 


( "l 

-2 / Bocosx\ 

Vio GRzJ 

^ 0.1 /iG J 


(7) 


Thus, if the line-of-sight component of the uniform magnetic field has a strength of ^ 0.1 /rG and we can 
observe at low frequencies such as 10 GHz, the magnitude of the rotation measure angle could be of the 
order of 10“^rad (« 6 °) for such a massive cluster. The observed scatter of the normalized rotation measure 


(RM = A( 75 /AQ[m]) is RM = 0 — 250 rad m ^ at the radius in the range of 0 , 


Ih ^ Mpc from the 


cluster center (Kim et al. 1991; Clarke et al. 2001). The value estimated in equation ([^ corresponds to 
RM 260 rad m“^ at the cluster center for the typical case of BqCOSx = 0.1 /xG. Our model thus seems 
consistent with the observations. This can be more quantitatively verified as follows. Using our model of 
the rotation measure angle (^) and the Press-Schechter mass function dn/dM ( Press, fc Schechter 1974 ) (see 
equation ( 0 )), we can compute the variance of the rotation measure for present-day clusters against the 
radius from the cluster center: 


(RM2)(r) = 




/M„ 


dM 
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A(p{r, z = 0) 
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( 8 ) 


In the calculation of the equation above, we use (cos^ x) = 1/3 for the ensemble average of angle x between 
the magnetic field direction and the line of sight, since the field orientation is considered to be random for 

each cluster. The solid line in Figure ^ shows the average profile ^ (RM^)(r) for Bq = 0.2/iG in the mass 
range of M > 5 x IO^^Mq. Note that this mass range corresponds to Lx ^ lO^'^erg s“^ for the X-ray 
luminosity in the ROSAT 0.1 — 2.4 keV band, where only the thermal bremsstrahlung emission is considered. 
Thus, this X-ray luminosity range is roughly consistent with the ROSAT sample of the Faraday rotation 
measurements used by Clarke et al. (2001). For comparison, square symbols represent the observational 
data from Clarke et al. (2001) with Tier errors arising from the Faraday contributions due to our Galaxy. 
It is clear that our model with Bq = 0.2/iG can roughly reproduce the observed magnitude and dispersion 
of the rotation measures as a function of the radius. In this sense, it is likely that we can at least estimate 
a correct magnitude of the Faraday rotation effect on the angular power spectra of CMB polarization fields, 
even though the uniform magnetic field is an unrealistic model, because the contribution to the angular 
power spectra comes from the second moments of the rotation measure angle by definition (see equation 
(I))- Furthermore, there could be an uncertain selection effect in the observational data. For example, one 
cluster sample with the large RM value of |RM| = 229 at the radius of r = 3 kpc has the X-ray luminosity 
of Lx = 2.0 X 10^2 erg s“^ or the mass of M « 7.0 x IO^^Mq, which is smaller than the lowest mass of 
M = 5 X IO^'^Mq used for the theoretical prediction shown by the solid line. To illustrate the sensitivity of 
the lowest mass to the theoretical prediction, the dashed line in Figure]^ shows the result of equation (||) for 
M > 10^'^Mq {Lx ~ 4.4 x 10^^ erg s“^) and Bq = 0.2/iG, and it clearly underestimates the observational 
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Fig. 1.— Angular profile of rms of the rotation measure for a cluster with M > 5 x as a function 

of the radius from the cluster center (solid line), which is computed using equation (|^). We hare considered 
Bq = 0.2/r G for the intracluster magnetic field strength. In order to demonstrate the dependence of the 
minimum mass used in the integration of equation (^, the result for M > is also shown by dashed 

line. The mass range of M > 5 x corresponds to the X-ray luminosity of Lx ~ lO^'^'^erg s“^ in 

the 0.1 — 2.4 keV energy band, which roughly corresponds to the ROSAT sample of the Faraday rotation 
measurements used by Clarke et al. (2001). For comparison, square symbols with error bars represent the 
observational data from Clarke et al. (2001), where the error arises from the Faraday rotation contribution 
due to our Galaxy. 


result. For the similar reason, it is thought that Clarke et al. (2001) derived a higher value of ^ 0.5 /rC for 
the uniform magnetic field model from their observational data. In addition to these facts, there is possible 
contamination of line emissions from the intracluster metal to the observed X-ray luminosity. Therefore, 
our model of Bq = 0.2 /iG often used in the following discussion could lead to the minimum magnitude of 
the Faraday rotation effect on the CMB polarization power spectra for a realistic effect. 


As will be shown in detail, the Faraday rotation effect generates secondary patterns of the CMB po¬ 
larizations expressed as Q ^ AipU and U ^ AipQ in terms of the Stokes parameters Q and U, where Q 
and U are those parameters of the primary polarization generated via Thomson scattering at the decoupling 
{z^ ~ 10^). Hence, equation (^ implies that, if we observe the CMB polarization through a cluster at longer 
wavelengths such as Aq ^ 6 cm (i^o ~ 5 GHz), the secondary signals could be comparable to the primary 
signals for massive clusters with BqCOSx = O.l/rC. In this case, we could directly measure the Faraday 
rotation effect in an individual cluster like the SZ effect on the CMB temperature fluctuations. This would 
be also true for particular clusters that have an unusual stronger magnetic field such as ^ 1 fj,G. Moreover, 
since the primary CMB polarization field has very smooth structure on angular scales of d 10' because 
of the Silk damping (^ilk 1968; and also see Hu, & Sugiyama 1995 in more detail), the angular profile of 










the rotation measure could be directly reconstructed in principle from the secondary signals generated on 
small scales even for a single observed frequency. These interesting possibilities are now being investigated 
in detail, and the results will be presented elsewhere (Ohno et al. 2001 in preparation). 


2.3. Angular po-wer spectrum of the rotation measure 

We are interested in the secondary effect on the angular power spectra of the CMB polarization due to 
the Faraday rotation effect caused by clusters. We thus need to compute the angular power spectrum of the 
rotation measure angle A(/j(0 ), which gives the strength of in the average sense for the all-sky survey at 
two-point statistics level. The angular Fourier transformation of the rotation measure ( 0 ) is given by 


Aipi{z) = 2Tr / dO0Aip{9,z)Jo{l0). 


(9) 


where 0vir = Vvir/dAiz), Jo{x) is the zero-th order Bessel function, and we have used the small-angle approx¬ 
imation (Bond & Efstathiou 1987). Since clusters are discrete sources, contributions to the power spectrum 
of the rotation measure are divided into the Poissonian contribution, C) , and the contribution from the 
spatial correlation between clusters, (Peebles 1980). From equation (^) the correlation contribution 

Cl depends on the ensemble average (cosyi cosx 2 )i where xi and X 2 are the relative angles between the 
magnetic fields and the line of sights for those two clusters, respectively, while depends on (cos^y), 
leading to the non-vanishing contribution for any distribution of y. If the origin of the magnetic field is 
primordial, the correlation between two intracluster magnetic fields could be non-vanishing even on separa¬ 
tion scales between clusters like that between the density fluctuation fields, and in this case we have some 
contribution from C) because (cosyiCOsy 2 ) ^ 0. However, it is known that the Poissonian contribution 
always dominates the correlation one because clusters are rare objects (Komatsu, &: Kitayama 1999|). For 


this reason, we here assume (cos yi cos y 2 ) = 0 and consider only the Poisson contribution C) ' for simplicity. 
Following the method developed by Cole & Kaiser (1988) and Limber’s equation, we can obtain an integral 
expression of 


C, 


.RM(P) 


n'‘°dz — 

Jo dz . 




M„, 


dM 


dn{M, z) 
dM 


\Aipi{M,z)(^ 


( 10 ) 


where V{z) and r{z) are the comoving volume and the comoving distance, respectively, and dn/M{M,z) is 
the comoving number density for halos with mass in the range of M and M + dM at redshift z. We compute 
dn/dM using the Press-Schechter model ([Press, fc Schechter 1974), and in this case we have 


dn{M, z) 

[2 p 

din a 

dM 

^ TT M 

dM 


Vn.e 


-CJ2 


( 11 ) 


where Vc = dc{z)/a{M) denotes the threshold, a{M) is the present-day rms mass fluctuation within the 
top-hat filter corresponding to the mass scale of M, 5c{z) « 1.69/Zl(z) is the threshold overdensity of the 
spherical collapse model, and D{z) is the linear growth factor of density fluctuations. 

Figure ^ shows the angular power spectra of the rotation measure through clusters for Bq = 0.2 
(solid line) and 0.5 fiG (dashed) assuming the observed frequency of vq = 10 GHz in the ACDM model. 
Note that we have again used (cos^ y) = 1/3. We have confirmed that dominant contribution to the power 
spectrum comes from the massive clusters with mass larger than lO^^M© at low redshifts such as z 0.5, 


similarly as the results of SZ effect due to clusters (e.g., Gooray 200C). 














- 9 - 






Fig. 2.— The angular power spectra of the rotation measure angle due the intracluster magnetic field 
calculated by equation & for the observed frequency of lOGHz in the LCDM model (flmo = 0.3, = 0.7, 

h = 0.7, riho = 0.05, (78 = !)■ We here considered two cases of the uniform magnetic fields Bq of 0.2 fj,G 
(solid line) and 0.5 fiG (dashed) strength across a cluster. Note that the power spectrum has the dependence 
on the magnetic field strength and the observed frequency as oc 


3. Secondary EfTect of Faraday Rotation in a Cluster on CMB Polarization 


In this section we develop a formalism for calculating the secondary Faraday rotation effect induced by 
magnetized intracluster hot plasma on the CMB polarization fields and their angular power spectra. 

The properties of CMB polarization fields can be described in terms of the Stokes parameters Q and U 
( Kosowsky 199(]| ; Kamionkowski, Kosowsky, & Stebbins 1997a,b; Zaldarriaga, & Seljak 1997). Let us consider 
quasi-monochromatic waves propagating in the z-direction, in which the electric and magnetic fields vibrate 
on the x-y plane. We can then express the electric field E in terms of the complex orthogonal unit vectors 
of circularly polarized waves as E = + £'Le*'^'^eL)e“*“‘, where = (e^ T *ey)/'s/2 and cr and 

Cl are the unit vectors of right-hand and left-hand polarized waves, respectively. Using those quantities, the 
Stokes parameters Q and U measured in the fixed (a:, ?/)-coordinates can be expressed as 


Q = \E-e,\^-\E-ey\^ = 2E^Ei^cos{5^-5i:), 

U = 2Re[{E-e,y{E-ey)]=2EnEi^sin{6K-5i^). (12) 


If there is an intervening magnetic field in the ionized plasma between the last scattering surface and the 
present epoch, the right-hand and left-hand polarized waves travel with different phase velocities by the 
Faraday rotation effect. Therefore, from equation we can derive the following equation to govern the 
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time evolution of Q and U in the presence of the magnetized plasma; 


dQ 

dt 

dU 

dt 


jj djSR, - ( 5 l ) 
dt 


—2ujbU, 


2ujbQ, 


(13) 


where 2ujb = d{SR — 5i,)/dt, t is the cosmic time and we have ignored other effects associated with the 
Thomson scattering. Note that ws is also related to the rotation measure angle /S.ip defined by equation i) 
as lob = d{Aip)/dt. Equation ( p^ clearly shows that the Faraday rotation causes the mixing between Q- and 
17-type polarizations. More importantly, even if the primary CMB anisotropies generate only the Q-type 
polarization in the statistical sense ( Zaldarriaga, fc Seljak 19^ ), the Faraday rotation mixing induces the 
[/-type polarization on the observed sky from the primary Q mode. 


Since we are interested in the secondary effect on the CMB polarizations, the solutions of equation ( 1 ) 
for Q and U along the line of sight with direction of 0 can be approximately solved up to the first order of 
as 


Q{to,9) s 

s Q{U,9)-2U{U,9) 

! dtujB{t,9) 

u{to,e) s 

s UiU,9) + 2Q{U,9) 

rto 

1 dtuiB{t,9) 


(14) 


where U and Q denote the primary U- and Q-polarization fluctuations generated at the decoupling epoch 
(z* « 10^), respectively, and we have neglected other secondary polarization fields generated at low red- 
shifts that have insignificant power compared to the primary signal on relevant angular scales for a realistic 
reionization history (e.g., Hu 2000). The quantities t* and to are the cosmic time at the decoupling and the 
present time, respectively. In what follows quantities with tilde symbol denote quantities at t*. 


The Stokes parameters Q and U depend on the choice of coordinates and are not invariant under the 
coordinate rotation. For this reason, it is useful to introduce the new bases of rotationally invariant polariza¬ 


tion fields, called the parity-even electric mode (E) and the parity-odd magnetic mode (B) (Zaldarriaga, & 


Seljak 1997; Kamionkowski, Kosowsky, & Stebbins 1997a,b; Hu, & White 1997a,b). Any polarization pattern 
on the sky can be decomposed into these modes. Since the secondary effect due to clusters of galaxies on 
the CMB polarization is important only on small angular scales such as 0 1°, we can safely employ the 


small-angle approximation ( Zaldarriaga, fc Seljak 1998 ). Based on these considerations, Q and U fields can 
be expressed using the two-dimensional Fourier transformation in terms of the electric and magnetic modes 
as 


^ [^[Cos20;-H;sin20;] e 

^ /(^ 


le 

•> 

le 


(15) 


where ipi is defined by lx + ily = Zexp(i^^), and Ej^ and are the two-dimensional Fourier coefficients for 
the E- and H-modes, respectively. The Fourier components Ei and Bi for the primary CMB polarization 
fields satisfy 


(Ei^,) = {2nfCj,i5\l-l'\ 

{^BiBl) = {2^fC^,5\l-l% (16) 
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where and are the angular power spectra of primary E- and B-modes, respectively. The primary 
B-mode polarization is generated by the quadrupole temperature anisotropies associated with the tensor 
(gravitational wave) and vector perturbations. If we assume the standard inflation-motivated model with 
a spectral index of rig = 1 for the scalar perturbations, the tensor and vector perturbations are sufficiently 
smaller than the scalar perturbations. In the following discussion we ignore the primary B-mode power 
spectrum for simplicity; = 0. 

From equation (|^ the two-point correlation functions of Q and U fields including the Faraday rotation 
mixing effect can be computed as 

Cq{ 0) = {Q{o)Q{e))^CQ{0) + 4Cfji0)c^^{e), 

Cui0) = {U{O)U{e))^Co{0)+4CQ{0)C^^{0), (17) 

with 

CQi0) = f-^^{Mio) + Mi0)), 

CuiO) = (18) 

where we have employed a reasonable assumption that the primary fields (Q and U) on the last scattering 
surface and the rotation measure held ub induced by low-z clusters are statistically uncorrelated. Note that 
(^RM(^) jg two-point correlation function of the rotation measure angle and can be expressed in terms 
of the angular power spectrum of the rotation measure as 

C^^{0) = (^Jdtu;Bit,0) jdt'ujB{t',e)^ = (19) 

where jg given by equation (|ic|). 


We are now in the position to derive the angular power spectra of E and B modes in the observed CMB 
sky. Following the method developed by Zaldarriaga, & Seljak (1998), Cei and Cbi can be expressed as 


cei = 7T [ 0d0[{CQi0)+cumJom + {CQ{0)-cumMio)], 
Jo 

Cbi = TT [ 0d0[{CQi0)+CumJom-iCQ{0)-CumJ4ilO)]. 

Jo 

Therefore, from equation (p^) we have 


Cei = C,,,+ Il'dl'W.fC 


'El 


ll'^EV 


Cbi = jl'dl'wS.Cj,,, 

where Wew and Wbw are the window functions defined by 




= 


2 / 0d0 \Jq{10)Jo( 1'0) - Ji{l0)Ji(l'0)] C^^{0), 
Jo 

2 [ 0d0[Jo{W)Jo{l'0) + J4{l0)J4{l'0)]C™{0). 


( 20 ) 


( 21 ) 


( 22 ) 


Equations (21) clearly show that the Faraday rotation mixing generates the non-vanishing angular power 
spectrum of the B-mode polarization from the primary E-mode. 
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4. Results and Estimations of the Detectability 


Throughout this paper, we employ the ACDM model with = 0-3, ^\o = 0-7, ^bo = 0.05, h = 0.7 
and CTg = 1.0 as supported from observations of the CMB anisotropies and the large-scale structure. As for 
the power spectrum of dark matter used in the calculation of the Press-Schechter mass function, we adopted 
the Harrison-Zel’dovich spectrum and the BBKS transfer function ( Bardeen et al. 1986 ) with the shape 
parameter from Sugiyama (1995). 


Figure ^ shows the window functions of B- and E-modes per logarithmic interval in I, and PW^,, 

for a given mode I' = 2000 as derived by equation (|^). We consider the magnetic field strength Bq = 0.2 
and the observed frequency 10 GHz in the ACDM model. One can see that the window functions for the 
mixing between E and B modes are well localized in I space. 


In Figure ^ we show the angular power spectra of E- and E-mode polarizations on the observed sky 
including the Faraday rotation mixing effect due to clusters for the observed frequency of lOGHz. The dot- 
dashed line shows the E-mode spectrum for Bq = 0.5/iG that corresponds to the value estimated from the 
observational data for a model of uniform intracluster magnetic field (Clarke et al. 2001), while the dotted line 
shows the result of Bq = 0.2/rG. For comparison, the E-mode power spectra with and without the secondary 
effect of Bq = 0.2/xG are also shown by solid and dashed lines, respectively. One can readily see that the 
generated E-mode power spectrum has a peak with magnitude of ~ 0.1 /rE(Eo/0.1 hG){i'q/10 GHz)“^ at 
I 1000. On smaller scales of I ^ 6000, the spectra of E and E modes coincide, and this is similar as the 
secondary effects caused by the peculiar velocity field of the ionized medium in the reionized universe (Hu 
200C; Liu et al. 2001 ). These results imply that, if the uniform magnetic field has a strength of ^ 0.5 ^G 
and we can observe at lower frequencies such as i/q iS 5 GHz, the magnitude of E-mode could be comparable 
with the primary polarization amplitude of ~ /iK. 


The detectability of the induced E mode polarization by the sensitive satellite mission Planck can be 
estimated following the method developed by Zaldarriaga & Seljak (1998). The relative error on the overall 
amplitude of the induced E component [3 can be estimated as 


f3 

where /sky is the fraction of observed sky and we here fix /sky = 0.8. As for the detector noise and the 
beam width, we assume w~^ = 4.72 x 10“^® and Og = Oiwhm/VS In 2 = 3.71 x 10“^ for the 30GHz channel on 
Planck, where 0fwhm is the full-width at half-maximum angle. Note that the 30GHz channel we here consider 
is the lowest frequency channel on Planck. From the dependence Aq of the rotation measure, the magnitude 
of the E-mode angular power spectrum at 30GHz becomes smaller by ^ 10“^ than the results estimated in 
Figure ^ Accordingly, for the magnetic held of Bq = 0.2 /iG, we have PjAp « 0.14. This means that this 
effect cannot be detected by the Planck mission for Bq = 0.2fiG. Gonversely, to get (3/A(3 1 at the 30GHz 

channel, the uniform magnetic held component must have ^1.5 ^G strength, and this seems an unrealistic 
value. However, we stress that even a null detection of this effect can set a new constraint on the magnetic 
held strength. It is also worth noting that this signal would be detected more effectively by a ground based 
experiments at lower frequency bands such as 10 GHz observing a small patch of the sky for a sufficiently 
long integration time to reduce the noise per pixel. 


/sky X]; (2^ + 1)/(1 ■ 


—1 

■ w '-e « 


ICbiY' 


(23) 
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Fig. 3.— This figure shows the weighted window functions l^W^, and defined by equation (p^ ) as a 

function of I for a given mode V = 2000. We here present their power per logarithmic interval in 1. We used 
the model of the magnetic field strength 0.2 p,G and the observed frequency 10 GHz. 



5. Discussions and Conclusions 

In this paper, we have investigated the secondary effect on the CMB polarization fields induced by the 
Faraday rotation effect of the magnetic field in an intracluster hot plasma at a low redshift. To illustrate the 
impact of this effect in the simplest way, we employed the simple model that the magnetic field has a uniform 
field strength such as 0.2 /iG across a cluster universally when the cluster formed, which is consistent with 
observations of the rotation measures in clusters (Clarke et al. 2001). As shown in Figure |^, our model 
can roughly reproduce the observed scatter of the rotation measures as a function of the distance from the 
cluster center (Kim et al. 1991; Clarke et al. 2001). In this sense, since the Faraday rotation effect on the 
angular power spectra of CMB polarizations comes from the second moments of the rotation measure angle 
(see equation (ill)), it is likely that our model can at least estimate a correct magnitude of the effect, even 
though detailed shapes of those spectra could be different for a realistic intracluster magnetic field as discussed 
below. We showed that the parity-odd H-mode angular power spectrum is generated on the observed sky 
by the Faraday rotation mixing effect as a new qualitative feature, when if the primary CMB anisotropies 
includes the parity-even A-mode only such as predicted by the standard inflation-motivated scenarios with 
pure primordial scalar perturbations. We estimated that the generated B-mode power spectrum has a peak 
with amplitude of ^JWCibJ^ 0.1 /rK(Bo/0.1 /rG)(i^o/10 GHz)“^ at I « 1000 under the plausible scenario 
of the cluster formation in the currently favored ACDM model. It was also shown that the lowest frequency 
30GHz channel of Planck can be used to set an upper limit of Bq ^ 1.5^G for a uniform component of 
the intracluster magnetic field. Detection or even null detection of the predicted B-type polarization will 
therefore be a new empirical tool to provide a calibration of an uncertain magnetic field in a cluster, combined 
with measurements of the gas distribution from the X-ray or SZ datasince. 

It is known that there are other secondary sources on the CMB polarization fields caused in the low 
redshift universe, and these nonlinear effects generally induce the B-type polarization pattern from the 
coupling with the primary B-mode. The main source is the gravitational leasing effect of the large-scale 
structure, leading to the B-type polarization that has a peak of ^ 0.3 /iK at I « 1000 and then a power of 
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Fig. 4.— The angular power spectra of E- and S-mode polarizations. The dashed and dotted lines show 
the E- and S-mode spectra including the Faraday mixing effect due to clusters at the observed frequency 
of 10 GHz for the magnetic field model of Bq = 0.2 /iG, respectively, while the dot-dashed line shows the 
H-mode spectrum for the case of Bq = 0.5 ^G that corresponds to the value estimated by Glarke et al. 
(2001) for a model of the uniform magnetic field. For comparison, the solid line shows the primary E-iaode 
power spectrum in the AGDM model. 


^ 0.01 fiK at I ^ 5000 for the similar AGDM model as considered in this paper ( Zaldarriaga, fc Seljak 1998 ). 
One of authors (N.S.) has quantitatvely shown that the secondary effect due to the peculiar motion of ionized 
medium in the large-scale structure also generates the H-mode polarization with ~ 0.01 fiK at 100 ^ I ^ 10"*’ 


for a realistic patchy reionization model of the universe (Liu et al. 2001). From these results, if the uniform 
magnetic field has a larger strength of i? ~ 0.5 /iG as a common feature of clusters and we observe at 
lower frequencies such as r'o ^ 5 GHz, the generated H-mode amplitude due to the Faraday rotation effect 
would be comparable with the primary polarization amplitude and larger than the other secondary signals 
at I ^ 1000, although observations at such low frequencies would suffer from large foreground contamination 


of the syncrhrotron radiation from our galaxy (Tegmark et al. 2000). Even if the uniform magnetic field 
is not as strong as 0.5 ^G, the sensitive multi low-frequency measurements will allow one to separate this 
signal from primary and those other secondary signals in principle thanks to the frequency dependence of 


cx i'q "* on Cbi, which is analogous to measurements of secondary temperature fluctuations by the SZ effect. 
There are possible contributions to affect our results that we have ignored in this paper. The intracluster 


tangled magnetic fields, as implied by the rotation measures of some sources inside the Coma cluster (Kim et 


al. 1990), would provide contributions to our results in addition to the contribution due to the uniform field. 
However, the coherent scale of tangled fields in a cluster is still unknown. The recent complete measurements 
of the rotation measures (Clarke et al. 2001) imply the 1 — 10 /rG strength of the tangled field under the 
assumption of the coherent scale of 10 — 100 kpc. It is also known that the strong radio halos embedded 
within a (cooling flow) cluster are pervaded by stronger magnetic fields on smaller scales such as ~ 50 fjiG 
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for Hydra A (Taylor, & Perley 1993; Taylor et al. 2001). Although the simplest tangled-cell model with a 
constant coherence length I has often been considered in the literature, the possible important mechanism 


to affect our results is the deporlarization effect ( Tribble 1991 ), leading to the cancellations to some extent 
between the Faraday rotations in independent cells with random orientations of the magnetic fields within a 
cluster. The finite beam size will also lead to the artificial cancellations between the Faraday rotations of cells 
covered within one beam. From these considerations, the depolarization effect may lead to an intuitive result 
that a feature of the H-mode spectrum can be sensitive to the characteristic angular scale, Ic, originating 
from the coherent length. This investigation is now in progress, and will be presented elsewhere (Ohno et al. 
2001). Recently, it has been suggested that outflow or jet acitivities of quasars or black holes leave behind 
an expanding magnetized bubble in the itergalactic medium (IGM) resulting in the possibility that the IGM 


is filled by magnetic fields to some extent ( Furlanetto, fc Loeb 2001 ; Kronberg et al. 2001 ). Even thgouh 
the IGM field has a smaller strength than ~ 10“® G on Mpc scales, a large volume coverage of the ionized 
plasma might lead to non-negligible contribution connected to the reionization history of the universe. 

Finally, we comment on implications of our results to the nonthermal processes in a cluster. The 
observations of diffuse synchrotron halos and nonthermal hard X-ray emission have implied the existence 
of nonthermal relativistic electrons in a cluster (e.g., Rephaeli, Gruber, fc Blanco 1999| ). These electrons 
are produced by the shock acceleration in an intracluster medium, and it will allow us to derive additional 
information on the physical conditions of the intracluster medium environment, which cannot be obtained 
from the thermal plasma emission only. The intracluster magnetic field should then play an important 
role to these nonthermal processes. In particular, it is being recognized that the nonthermal process can 
be a unique probe of measuring the dynamically forming clusters before thermalization of the intracluster 
gas. Recently, Waxman & Loeb (2000) pointed out that the synchrotron radiation from those forming 
clusters could be foreground sources of the CMB temperature fluctuations at low frequency observations 
such as 10 GHz. Those clusters could also produce the high energy gamma ray emission through the inverse 
Compton scattering of CMB photons by relativistic electrons and thus be identified as gamma ray clusters 
by the future sensitive telescope GLAST, which are difficult to detect by X-ray or optical surveys because 
of their extended angular size of about ^ 1° ( Totani, fc Kitayama 2000 ). Based on these ideas, it will 
be very interesting to investigate cross-correlations between the synchrotron radiation or the gamma ray 
sources and the Faraday rotation signal in CMB polarization map. This investigation is expected to provide 
complementary constraints on the nonthermal physical processes in the forming clusters. 
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